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ABSTRACT 



T. Belloni 1 



We report on the results obtained for the millisecond transient X-ray pulsar XTE J1807-294 in a 40 
min orbital period system, based on an XMM-Newton ToO observation carried out during March 2003. 
The source was found at a luminosity level of about 2 x 10 36 erg s _1 in the 0.5-10 keV range (assuming 
a distance of 8 kpc). We confirm the presence of the 5.25 ms pulsations (after accounting for the orbital 
modulation) and find a pulsed fraction of 5.8% in the 0.3-10 keV band. The pulse shape in nearly 
sinusoidal. The spectral continuum of the source is well fitted by an absorbed Comptonization model 
plus a soft component. No emission or absorption lines have been detected in the 0.5-10 keV range with 
upper limits of 10-40 eV. The reported analysis represent the first detailed study of this source, the 
fourth belonging to the ultra-compact binary system class hosting an accreting neutron star. 

Subject headings: accretion, accretion disks 
stars: neutron 



binaries: close — star: individual (XTE 1807-294) 



1. INTRODUCTION 



In the past few years conclusive evidence that Low 
Mass X-ray Binaries (LMXRBs) contain fast spinning, 
weakly magnetic neutron stars has been gathered. Coher- 
ent periodicities have been discovered in four sources SAX 
J1808.4-3654 (401 Hz, Wijnands & van der Klis 1998), 
XTE J1751-305 (435 Hz, Markwardt et al. 2002), XTE 
J0929-314 (185 Hz, Galloway et al. 2002) and recently 
XTE J1807-294 (191 Hz, Markwardt, Smith & Swank 
2003). 

All these four sources are transients, i.e. exhibit only 
sporadic activity and for most of the time remain in a 
state of low level activity (for a review see e.g. Campana 
et al. 1998a). These transient pulsars share many similar- 
ities among them, but are rather peculiar among transient 
sources containing a neutron star (usually named Soft X— 
ray Transients, SXRT). Their peak X-ray luminosities are 
rather faint (~ 0.01 L^dd), whereas SXRTs have peak lu- 
minosities of 0.1 — liEdd- They show X-ray pulsations, 
whereas for SXRTs and LMXRBs tight upper limits ex- 
ist (e.g. Vaughan et al. 1994). The discovery of X-ray 
pulsation during the outburst phase has also allowed to 
reveal an orbital periodicity, which in all cases is shorted 
than - 2 hr (2.01 hr for SAX J1808.4-3654, Chakrabarty 
& Morgan 1998; 42.4 min for XTE J1751-305 Markwardt 
et al. 2002; 43.6 min for XTE J0929-314 Galloway et al. 
2002; 40.1 min for XTE J1807-294, Markwardt, Juda & 
Swank 2003). These periods are much shorter than what 
usually found in SXRTs ( 4 hr) . These differences led 
some authors to suggest that transient pulsars might form 
a separate (sub-)class among SXRTs, called faint tran- 
sients (Heise et al. 1998; King 2000; in't Zand 2001). 

Faint transients are poorly known. The best studied 
source is SAX J1808. 4-3658. It was studied across three 
outbursts (Wijnands & van der Klis 1998; Wijnands et al. 
2001), in quiescence (Campana et al. 2002), in the opti- 
cal (Homer et al. 2001) and radio (Gaensler et al. 1999) 
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bands. 

Concerning XTE J1807-294 very little is known. The 
source was discovered during a Galactic plane scan by 
RXTE/PCA (Markwardt, Smith & Swank 2003). The 
source was observed at a level of 33, 38, 58, 41 and 20 
mCrab (2-10 keV) on 2003 Feb 16.7, 19.8, 21.6, 22.6 and 
Mar 13.9, respectively. A follow-up Chandra observation 
on 2003 Mar 10.8 allowed to obtain a much better position 
(1" uncertainty). Pulsations were still detectable and al- 
lowed a refined orbital solution to be obtained from RXTE 
data (Markwardt, Juda & Swank 2003). 

Here we report on an XMM-Newton (discretionary 
time) TOO observation of XTE J1807-294 on 2003 Mar 
22 (see also Campana et al. 2003). In Section 2 we de- 
scribe the observational data and the reduction strategy. 
In Section 3 we describe our results on spectral and timing 
analysis. Conclusions are reported in Section 4. 

2. XMM-NEWTON OBSERVATION AND DATA ANALYSIS 

The XMM-Newton Observatory (Jansen et al. 2001) 
includes three 1500 cm 2 X-ray telescopes each with a Eu- 
ropean Photon Imaging Camera (EPIC, 0.1-15 keV) at 
the focus. Two of the EPIC imaging spectrometers use 
MOS CCDs (Turner et al. 2001) and one uses pn CCDs 
(Striider et al. 2001). Reflection Grating Spectrometers 
(RGS, 0.35-2.5 keV, den Herder et al. 2001) are located 
behind two of the telescopes. 

XTE J1807-294 was observed on March 22, 2003 from 
13:57 to 18:40. MOS1 was operating in Small Window 
mode, MOS2 in Prime Full Window mode and pn in tim- 
ing mode, all with thick filters. RGS1 and RGS2 were 
operating in the Spectro + Q standard mode. Due to a 
problem in the field acquisition the source was originally 
pointed 1 arcmin apart. This problem was then corrected, 
resulting in an overall efficiency decrease of the pn, MOS1 
and RGS instruments. We generated final product using 
sas v. 5.4.1. A small soft proton flare event occurred 
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at the beginning of the observation, we excluded the first 
2.5 ks from the analysis. We obtained net exposure times 
of 14 ks and 9 ks for MOS2 and the other instruments, 
respectively. 

2.1. MOS cameras 

The source is extremely bright such that the MOS cam- 
eras are saturated and the radial support structure of the 
mirrors is visible. Piled-up rate for the MOS2 camera is 
~ 8 counts s _1 (0.5-10 keV). Source position is consistent 
with that determined with Chandra (Markwardt, Juda & 
Swank 2003). 

A simple way round to pile-up is to extract photons 
for spectral analysis in an annulus around the source, thus 
excluding the inner (piled-up) core. We followed this strat- 
egy for analysing the MOS2 data, whereas for MOS1 data 
the window is too small (see also the analysis carried out 
on Mkn 421, presenting similar problems, by Ravasio et 
al. 2003, in preparation). We extracted photons from an 
annulus of 30" and 60" inner and outer radii, respectively. 
Background was extracted from two different outer circu- 
lar regions 2' each. We selected only single-pixel events 
and make sure that no pile-up is present (through the sas 
tool epatplot). This is also in line with prescriptions by 
Molendi & Sembay (2003) 4 . The count rate within the 
0.5-10 keV energy band is 1.01 ± 0.01 counts s _1 . Re- 
sponse matrix and ancillary files were then generated with 
rmf gen and arf gen. Spectral data were rebinned to have 
80 counts per bin. 

2.2. pn camera 

Timing mode was selected to have the opportunity to 
reveal X-ray pulsations. Given the high count rate of the 
source, this mode is particularly useful to avoid any pile- 
up. We extracted the source spectrum from the raw image, 
taking a column 10- pixel wide and the background from 
two separate boxes 5-pixels wide. We consider only "sin- 
gle" and "double" events. No X-ray bursts were detected. 
The source count rate in the 0.5-10 keV band is 34.0 ±0.1 
counts s _1 . Background is at a level of 1.4% in the full 
band. Pile-up is negligible and dead time is at a level of 
~ 1.5%. Response matrix and ancillary file for spectral 
analysis were generated using sas tasks. Spectral data 
were rebinned to have 80 counts per bin. 

2.3. RGS 

The standard sas procedure rgsproc was used to de- 
rive the first order RGS spectra. In the 7-36 A band rates 
are 0.39 ± 0.01 and 0.56 ± 0.01 count s" 1 for RGS1 and 
RGS2, respectively. For spectral analysis together with 
the other instruments we rebinned the data to 100 counts 
per bin. We further inspect the RGS spectra for emission 
or absorption lines. In this case we rcbin the spectra to 20 
counts per bin (see below). 

3. SPECTRAL ANALYSIS 

The spectral analysis has been carried out in the 0.5- 
10 keV energy range for MOS2 and pn and in 7-36 A for 
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the two RGSs, using XSPEC (vll.2). All spectral un- 
certainties are given at 90% confidence level for one de- 
gree of freedom (A^ 2 = 2.71). We first tried the stan- 
dard model for SXRTs, i.e. an absorbed soft compo- 
nent (black body) plus a hard component (power law; 
e.g. Campana et al. 1998b). The fit is good with a 
Xred = 1-10 (1571 degrees of freedom, d.o.f.). Fit pa- 
rameters can be found in Table 1. The column density 
is (6.3 ± 0.1) x 10 21 cm~ 2 which is slightly higher than 
the Galactic value of 3 x 10 21 cm~ 2 . A comparable (but 
more physical) model is obtained by replacing the power 
law with a Comptonization model like C0MPTT (Titarchuk 
1994) as observed in SAX J1808.4-3658 (Gierlinski, Done 
& Barret 2002 but not in XTE 1751-305, Miller et al. 
2003). In this case we obtain x 2 oc i = 1-10 (1569 d.o.f., see 
Fig. 1). Systematic errors at 1.6% level provide an ac- 
ceptable fit. The column density is 4.6t°;2 x 10 21 cm~ 2 , 
consistent with the Galactic value. The black body tem- 
perature is k T — 0.80 ± 0.03 keV and its radius (assuming 
spherical symmetry) is R = (2.2 ± 0.2) d\ f 2 km (with dg 
the unknown distance to XTE J1807-294 in units of 8 kpc 5 
and / the spectral hardening factor obtained as the ratio 
of the color temperature to the effective temperature, e.g. 
Merloni, Fabian & Ross 2000). The temperature of the 
soft input photon is 0.24 ± 0.02 keV. In the Comptoniza- 
tion model the temperature of the corona and the plasma 
optical depth are tightly related. The best fit is for a rel- 
atively cool corona (kT = 9.6^° 2 keV, i.e. kT > 5.6 
keV) and an optical depth (r < 4.2). A very hot corona 
(kT > 370 keV) with a negligible optical depth (r ~ 0.02, 
r < 2.68 at 90% confidence level, c.l.) is also (formally) 
acceptable. However, with a detector sensitive up to 10 
keV, it is impossible to test the presence of such a high 
coronal temperature. The unabsorbed 0.5-10 keV flux for 
our best fit spectrum amounts to 2.6 x 10~ 10 erg s _1 cm~ 2 . 
Source luminosity is 2.0 x 10 36 c?| erg s _1 . The extrapo- 
lated luminosity to the 0.01-100 keV band is a factor ~ 1.7 
larger. 

The hard component contributes almost the entire (un- 
absorbed) flux (87%). The black body component is de- 
tectable only at low energies. Other soft components have 
been tried with similarly good results. We cannot discrim- 
inate between an emission from a black body or from an 
accretion disk (Td ~ 0.2 keV, Rd ~ 10 km). A neutron 
star atmosphere (Gansicke, Braje & Romani 2002) and 
a thermal spectrum however provided a worse fit to the 
data. 

3.1. Line search 

We searched for an iron lines in the pn and MOS data 
with negative results. Taking the line energy in 6.40-6.97 
keV and 0.1 keV (best fit) or null values for the line width, 
we derive an upper limit on the equivalent width of 18 and 
25 eV, respectively. 

Following Miller et al. (2002) we searched the RGS spec- 
tra for emission or absorption lines. We divided the spec- 
tra (rebinned at 20 counts per bin) in slices of 3 A and 
inspect them for lines in the energy band 7-20 A. We fit- 
ted the sliced spectra with an absorbed (VARTBABS, Wilms, 
Allen & McCray 2000) power law model. All fits are ac- 
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ceptable. No emission or absorption lines are found. We 
directly searched for Ne K edge line at 14.25 A following 
the evidence for enhanced absorption in ultra-compact bi- 
naries (Schulz et al. 2001; Juctt, Psaltis & Chakrabarty 
2001). No Ne edge is visible and the fit procedure is not 
sensible to its presence. We also searched for Ne Lya line 
at 12.13 A (in the range 11-13 A). No line is present and 
an upper limit of 41 eV can be set (95% c.l.). In case 
of a narrow line the limit is 9 eV. Detected Ne lines in 
4U 1626-67 had equivalent witdhs of 20 A (Schulz et al. 
2001). 

4. TIMING ANALYSIS 

We used the sas task barycen to correct the pn data. 
The time resolution is 29.56 ^s (Kuster et al. 2003). A 
significant signal in the power spectrum is readily visible 
in the data (Campana et al. 2003). However, the power 
spectrum shows a large number of peaks likely resulting 
from aliasing of the short orbital period (see Fig. 2). The 
correct pulse period can only be obtained after correcting 
for the orbital modulation. We assumed the orbital pe- 
riod derived by Markwardt, Juda & Swank (2003) using 
RXTE data, i.e. P or b = 40.0741 min, and searched for 
the projected orbital radius and time of 90° mean longi- 
tude assuming the hypothesis of circular orbit. Stepping 
on those two parameters and, for each case, folding the 
events into a 10-bin light curve. The presence of a modula- 
tion in the folded light curve was then estimated with a y 2 
test. This is similar to the analysis carried out by Kirsch & 
Kcndziorra (2003). We successfully recover the pulsation 
with a projected orbital radius a x sini = 4.7^qq 4 and an 
epoch of 90° mean longitude T = 52720.68644(13) (90% 
c.l., see Fig. 3). Pulsations were detected at 5.2459426(2) 
ms (90% c.l.). 

After correction for the orbital motion and in the full 
band used for timing analysis (0.3-10 keV) the (back- 
ground subtracted) pulsed fraction amounts to 5.8 ±0.3% 
(90% c.l., without correction the pulsed fraction is only 
1.9%). The fit is not completely satisfactory (Xred = 4-8). 
An additional sinusoidal component at half of the spin 
period is highly significant (F-test probability 3.5 a and 
X? od = 1.2). The pulsed fraction of XTE J1807-289 is 
comparable to what observed in SAX J1808.4-3658 (4%, 
Wijnands & van der Klis 1998), XTE 1751-305 (4%, Mark- 
wardt et al. 2002) and XTE J0929-314 (~ 5%, Galloway 
et al. 2002). The pulsed fraction increases with energy, 
going from 4.3 ± 0.7 in 0.3-2 keV to 7.4 ± 0.7 in 2-5 keV 
and to 7.5 ± 1.4 in 5-10 keV. Also for these curves two 
sinusoidal components provide a much better fit. 



The orbital modulation is not directly visible in the 
data. Folding the barycentered pn light curve at the 
orbital period reported by Markwardt, Juda & Swank 
(2003), we can put a 95% upper limit on the modulated 
component of 1%. 

5. CONCLUSIONS 

We report on XMM-Newton observation of the fourth 
millisecond low mass X-ray pulsar XTE J 1807-294. We 
reveal the source still in a bright state ~ 30 d after its 
discovery. The source luminosity is 2 x 10 36 erg s _1 at a 
fiducial distance of 8 kpc. 

Our results are very similar to what observed in XTE 
J1751-305 (Miller et al 2003). The hard spectral com- 
ponent contributes most of the observed flux (87%), even 
though a soft component (a black body) is needed by the 
data. The equivalent radius of the black body component 
is small, R = (2.2 ± 0.2) d\ f 2 km, suggesting an emitting 
hot spot onto the neutron star surface. 

No emission or absorption lines are observed in the data. 
We study in more detail Ne lines following the sugges- 
tion by Yungelson, Nelemans & van den Heuvel (2002) 
and Bildsten (2002) that the companions of faint tran- 
sient sources are Ne-rich white dwarves. We did not find 
any evidence of emission lines (41 eV and 9 eV 95% upper 
limits for a broad or narrow line) or absorption edges. Iron 
lines are not detected either with upper limits of ~ 25 eV. 

The hard component (the main contributor to the source 
flux) can be better described with a Comptonization model 
rather than a simple power law model. This is at variance 
with the case of XTE J 175 1-305. The optical depth and 
the plasma temperature are highly covariant. In any case a 
firm upper limit on the optical depth of 4.3 (90% c.l.) can 
be set. This is in line with the predictions of Titarchuk, 
Cui & Wood (2002). They suggested that X-ray pulsa- 
tions can be observed in sources with an optical depth 
lower than <~ 4 at a plasma temperature of ~ 20 keV (as 
in our case). For sources with a higher optical depth, e.g. 
t = 11.7 ± 0.4 in the Z source GX 349+2 (Di Salvo et 
al. 2001) or t = 5.3 ± 0.6 in the atoll source 4U 1728-34 
(Di Salvo et al. 2000), pulsations would be more easily 
suppressed by scattering. 

We are grateful to the XMM-Newton project scien- 
tist, F. Jansen, for granting time to this unsolicited tar- 
get of opportunity observation. We thank the referee for 
her/his helpful comments. This work is partially sup- 
ported through ASI and Co-fin grants. 
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Table 1 

Spectral fit with an absorbed black body component plus a power law or a COMPTT component. 



Parameter 



Value 



6.3±tj"i x 10 21 cm" 2 
0-81±g-gt keV 
l.8t° Q 2 2 km 
1 Q6+ 002 
1.10 (1569 d.o.f.) 



Black body T 
Black body R* 
Photon index 

Xrcd 



4.6±^ x 10 21 



Black body T 
Black body R* 
Input To 
Plasma T p 
Optical depth 

Xrcd 



cm 



0.80±°;°ij keV 



2.2 



+0.2 



km 



0.24t°°| keV 
.+102 keV 



9.6 



4.0 

1.10 (1569 d.o.f.) 



Note. — * Black body radius at a distance of 8 kpc. 
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Fig. 1.— XMM-Ncwton spectra of XTE J1807-294 with overlaid the beEltQisg-)fcd(rJfa@c\o^dy plus COMPTT model. From top to bottom 
there are the pn and the MOS2 spectra. In the lower panel are shown the residuals in terms of \ 2 - 
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Fig. 2. — The upper panel shows the contour plot in the a x sini - phase plane of the best correction for the orbital motion (assuming 
the orbital period of 40.0741 min, Markwardt, Juda & Swank 2003). Contours (from outside) refers to 99.73%, 95%, 90% and 68.3% c.l., 
respectively (the best fit is marked with a small line). Orbital phases are in arbitrary units (from the time of the first photon received 
by XMM-Newton) . In the lower panel the signal power (in terms of x 2 ) is shown. The abscissa is in units of /is from the best period 
Pbcst = 5.2459426 ms. 
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Fig. 3. — Pn light curve of XTE J1807-294 folded at the neutron star spin period. From top to bottom, these refer to 0.3-10 keV, 0.3-2 
keV, 2-5 keV and 5-10 keV, respectively. The curve is repeated twice for clarity. 



